Suppressing the large scale curvature perturbation by interacting fluids 
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The large-scale dynamics of a two-fluid system with a time dependent interaction is studied 
analytically and numerically. We show how a rapid transition can signihcantly suppress the large- 
scale curvature perturbation and present approximative formulae for estimating the effect. By 
comparing to numerical results, we study the applicability of the approximation and find good 
agreement with exact calculations. 
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I. INTRODUCTION 

There are many instances in cosmology, when the large-scale evolution of perturbations [l| in a multi-component 
system needs to be solved [2|. The use of curvature and entropy perturbations has become a de facto standard in 
[ — \ this formalism, pioneered by Kodama and Sasaki [3]. In a multi-component system, interactions between fluids are 
I an important aspect in determining the evolution of the curvature perturbation Q • The importance of such systems 
' is evident: interacting fluids are the cornerstone of many iinportant cosmological scenarios, such as reheating at the 
^ ; end of inflation @, 0,1, i, [O, El or the curvaton scenario [11 [H, [H, [H, |I3 ■ 

' The total curvature perturbation, C, generally evolves in multi-component fluid whenever the non-adiabatic pressure 
is non-zero. In particular, if the fluids are interacting the total curvature perturbation changes. Whether the evolution 
changes the amplitude of the large-scale curvature perturbation significantly or not, depends on the exact nature of the 
' system. Efficient amplification, or damping, of ^ on large scales by later evolution allows one to modify the previously 
produced perturbation spectrum, relaxing constraints on the scale at which the perturbations were produced. A 
natural framework for such mechanisms is e.g. within a traditional multiple inflationary scenario or a string 
, landscape picture [H, HI . 

The transfer of energy between fluids can be described by a number of ways: for example in Q a constant interaction 
T— I , between different components is assumed whereas authors of [111 and [13] use the so-called sudden decay approximation 
■ instead. In this approximation, the two fluids evolve independently until energy is transferred rapidly from one fluid 
to the other at a particular point in time. This method allows one to estimate the magnitude of the total curvature 
perturbation analytically. Here we relax these assumptions and allow for a time-dependent interaction while evolving 
the full large-scale perturbation equations with no further approximations. Research in this direction has recently been 
presented in [Tlj and [l3| , where locally different decay rates of the inflaton are generated by spatially varying reheating 
' temperature and couplings. Our approach is to study how turning on the interaction affects the perturbations and 
the significance of the transition time scale. When the interaction is turned on rapidly so that it is well modeled by 
O ■ a step-function, we can proceed with analytical methods and derive a useful approximation to the numerical results. 
^ ^ In the fluid approach with constant interaction term, the microphysics cannot be properly accounted for as the fluid 
^2 begins to decay simply when its decay width is of the order of the Hubble rate, V ^ H. Such a simple description does 
not always properly reflect the physical process, however, since there can be additional physical scales and processes 
involved that can play an important role. Cosmological processes involving different physical scales, such as phase 
transitions, do not need to conform to the simple fluid description but are better described by a time-dependent 
^ I , interaction. 

' In this paper we will consider time dependent interactions between different fluids. In section II we present the 
governing equations of perturbations in the Newtonian gauge. In the following section we compare the evolution of 
perturbations by means of numerical and analytical methods. We present analytical formulae which give an accurate 
approximation of the behavior of curvature perturbation variables and study the validity of these approximations. 
We end the paper with a short conclusion. Precise calculations of the evolution of the curvature perturbation of the 
decaying fluid during transition are presented in the Appendix. 
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II. PERTURBATION EQUATIONS 



Following the notation of [l|, 01 : we consider linear scalar perturbations about a spatially flat Fricdmann-Robertson- 
Walker -background in the Newtonian gauge: 

- ds^ = g/.^dz'^dx'' = (1 + 2(t>)dt^ - a(t)^(l - 2^)5,jdx^dx\ (1) 

where a is the scale factor. The background evolution is determined by the Einstein's equations, G^n = 87:GTfj_,y, and 
the continuity equations of the individual fluids: 

/jo(Q) +3iJ(l + W(„))po(Q) = Qq, (2) 

where is the equation of state and Qa describes the interactions between the different fluids. Total energy density 
is conserved, hence J2a ~ ^- Perturbing the covariant continuity equations, one finds the evolution equations of 
the perturbed energy and pressure densities (on large scales) 0|: 

SPa + 3H{Sp(^a) + SP(a)) " HPo{a) + ^O(a))^ = Qa(f> + SQa- (3) 

In addition, we have the Gq component of perturbed Einstein's equations 

3H{ijj + H(j)) = -AirGSp, (4) 

where 6p = J2a^Pa- Note that for perfect fluids (f) = ip in the Newtonian gauge, and hence given the equations of 
state, W(q,), and the interactions between the fluids, Qa, one can evolve the individual fluid perturbations along with 
the metric perturbation 0. 

The curvature perturbation of the fluid component a is defined to be the metric perturbation on uniform a-fluid 
density hypersurfaces, 

Cia) = -'ip-H- = -V' ; , (5) 

PO(.a) Po(a) 

where comma represents derivative with respect to the number of e-folds N, ' = d/dN = d/d(hia). The total 
curvature perturbation can be expressed as sum of the individual curvature perturbations, 

c-E^C(.), (6) 

a 

where po = J2aPo{a)- 

From the definition of the curvature perturbations (O one discerns that if Po(q) ~ 0> tlie variable becomes ill- 
defined. Fortunately, one can always use the equation for the total curvature perturbation, ^ and hence the problem 
can be circumvented if only one of the individual curvature perturbations becomes ill-defined. If more than one 
component suffers from the same problem, one can use an alternative formalism, e.g. by means of individual density 
perturbations Sp(^a)i or simply choose the set of evolving variables suitably before each problematic point. 

Here we consider an interacting two-fluid system, where energy density flows from one fluid, pi, to another, p2'. 
Qi = — r/(iV)po(i) , Q2 = r/(iV)po(i). During the evolution of the system, the curvature perturbation of the second 
fluid, C(2) becomes ill-defined at some point during the period of interest. Therefore, we choose our set of dynamical 
variables such that the system is described by the metric perturbation, (f>, the total curvature perturbation, C, and 
the curvature perturbation of fluid one, The equations of state of the fluids, are taken to be constants for 
simplicity. 

In terms of theses variables, the relevant background equations ^ take the form: 

P'o = -3(po + t^(i)Po(i) +t^(2)(Po-Po(i))) (7) 
P'on) = -(^3(l + ^(i)) + ^^)po(i), (8) 
The equations of motion of the perturbations can be derived from the definition of Q^j and Eqs. Q and ^ 

0' = ^{C + <f>)-4>, (9) 

2po 

C = -^(u;(2)-c.(i))(C-C(i)), (10) 
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The function f{N) G [0, 1] describes the (time-dependent) coupUng between the two fluid components. Physically, 
f{N) sets how quickly the interaction is turned on up to the full rate F. The form of the function f{N) has to be 
determined by microphysical properties of the interacting fluids together with necessary information on the evolution of 
the whole system. For example, the interaction strength can naturally be temperature dependent or phase transitions 
can significantly change the properties of the interacting system. To model and estimate the importance of such effects 
in general, we here adopt a simple Ansatz for the interaction without addressing the question of its origin. From the 
equations ([9|)- (fTT|) we see how a varying coupling, f'{N) ^ 0, between the fluids leads to an additional coupling of 
the metric perturbation, (j), to the curvature perturbations, leading to new effects. 



III. TIME-DEPENDENT INTERACTION 



As a first approximation, we choose the transition to be proportional to the Heaviside function, i.e. f{N) = 
9{N — iV*), with being the time when the transition begins and the strength of the interaction is denoted by 
F. This choice has the advantage that it can be solved analytically. The fluids are allowed to interact without any 
further approximations, in contrast to the sudden decay approximation, where the energy transfer between fluids 
occurs instantly when T/H reaches a critical value. 

With this choice of f(N), we can integrate the equation of motion of C(i), Eq. pip . The integration requires some 
care and the technical details are presented in the Appendix. The result is 

2 + x.ln[^] + ^'^^ 

where = (3H(A^*)(1 + cj(i)))/r and C(i)± = C(i)(-^*i) are the values right before and after the discontinuous jump. 

Physical processes however rarely happen as quickly as the Heaviside function would require. Therefore a steep 
hyperbolic tangent is a more realistic choice to model f{N). From the point of view of the dynamics of our system 
e.g. Eq. (|lip this means that during transition, when the term proportional to /'(N) dominates, the system is driven 
towards a equilibrium where Qi) = —0. Now the magnitude of F tells if the system will reach the equilibrium. We 
can now compare the magnitude of the jump from equation p2p to that of —0*, i.e. if Ci+ — Ci_ from eq. p2p is 
larger than —(f)^ our jump-formula is longer valid. This comparison leads to a simple equation for a;*, which can be 
solved numerically. The result is that eq. (fT^ can be used when F/i7(iV*) <C 11.76(1 -|- 'j-'(i))- When we increase F 
the system reaches the equilibrium and C(i)+ — —4>{N*) + Ci-- 

In order to determine the value of the metric perturbation cf) and the Hubble parameter H = at = 
(where SttG/S = 1) we need to evolve rest of the equations of motion imtil the beginning of the transition. Before 
the transition, Qj^) is clearly a constant and hence C(i)- = C(i)(-^o)- Therefore the equation of motion of C can be 
separated and integrated: 

/■(m (1+^(1))P0(1)(A^0) + (1+^(2))P0(2)(A^0) /^.,, . (N\^^r (N\ 

(1 + W(i))Po(i)(^o)e'' + (1 + W(2))Po(2)(Ao) V / 

where Nq denotes initial time, and we have defined AA^ = N — Nq and Auj = uj(^2) ~ ^(i)- 

Unfortunately, we are unable to solve for ipi^) in the general case and hence we have concentrated on the two 
special cases where one of the fluids is initially dominating the energy density of the universe. 

A. Dominant decaying fluid: po(i)(^*) ^ Po(2){N*) 
Before the transition different densities evolve according to 

Po(A^) -Po(i)(A^o)e-3(^+"<^')^^+ (po(A^o)-Po(i)(iVo))e-3(i+-(^))^^, ^^^^ 
Po(i)(A)=po(i)(A^o)e-^('+"<^')^^. 
Since /0(o)i dominates when N < N^, equation p3p implies that 

C(A^) = {({No) - C(i)(A^o)) e-3A-AAr ^ (^j.^ 
Inserting this into the equation of motion of (j) and solving the resulting first order differential equation, we have 
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9w(i) — 6lJ(2) + 5 V ^ ' / 5 + 307(1) 

+ kiVo) '^'T'^^L fCW - C(i)(^^o)) + ^^^^;^C(i)(^o)le-('^+-<^))^-/^. (16) 
L 9a;(i) — dw(2) + 5 V ^ ' / 5 + 3a;(i) ^ ' J 

Eqs. and (fT6|) along with the fact C(i)- = C(i) (^o) tell us the magnitude of the jump in the curvature perturbation 
of fluid one, C(i)+ ~C(i)-- These equations reveal that the magnitude of the jump depends crucially on the ratio F/i/*, 
in addition to the initial values. Note further that in most cases where a;(2) — (^(i) > or the initial state is adiabatic, 
(piN) reaches its asymptotic value, 3(1 + a;(i))C(i)(-/Vo)/(5 + 3w(i)), quite rapidly. Interestingly, non-adiabatic initial 
states with a;(2) — W(i) < may lead to exponentially amplified curvature perturbations. 

In addition to the magnitude of the jump in ((^i) , an interesting quantity is the final value of the total curvature 
perturbation, Ci^oo)- The magnitude of C also makes a (continuous) jump during the transition, but reaches a 
constant value much more quickly than This fact allows us to derive quite a useful approximation for the final 
value of the curvature perturbation. 

As the energy flows quite rapidly into Po(2)j C becomes a constant almost immediately after the transition. In 
addition, since Po(i) vanishes quickly, ^ changes more rapidly than and hence during this time we can approximate 
C(i) ~ C(i)+- Similarly, since H'{N)/H[N) = —3(1 + a;)/2, with uj = (aj(i)po(i) + ^{2)Po{2)) I Po being the effective 
equation of state, the Hubble parameter changes slowly compared to /9o(i) and hence H{N) ~ H{N^) = during 
the transition. 

Based on these arguments we can estimate the difference between the relative rates of change of Po(i) and p^: 

^ - - = -3(1 + ^(1)) - 77^ + 3(1 + c.) ^ --^ + 3Ac.„ (17) 
Po(i) Po H[N) 

where Aw, = lo^ — oj^i^^ and uj^ = lu{N^). Assuming that T/H^ > 3Aaj, and integrating both sides of equation (jl7p . 
we obtain 

Po(i)(iV) ^ exp[(3Ac., - ^){N - N^)]po{N). (18) 

Using this in the equation of motion of (, Eq. (fTU|) . we get 

C ^^oo = C(i)+ - (C(i)+ - C* 7T— ^^P'^TTi r ' 19 



where = C(-^*)> Pi* ~ Po(i)(-^*) and — po(A^*)- When the system in initially dominated by decaying fluid, 
we may set ~ cji. This equation along with Eqs. (jl2p and (jl3p determine the flnal value of the total curvature 
perturbation in terms of initial conditions and interaction strength ratio T/H^,. 

B. Sub-dominant decaying fluid: p(o)i(A'^*) ^ p(o)2(A^*) 
Supposing that the second fluid dominates before transition, i.e. Po(i) Po(2) we have 

poiN) ~ po(2)(iV) = Po(2)(iVo)e-=^(i+"<^')'^^. (20) 
Hence it follows from equation ([12]), that 

ciN) = am), (21) 

regardless of the initial conditions of C and The corresponding solution of (f> is 

C^iN) ^ UiNo) + 3(1 + -(^))C(^0)\ -(3.<.>+5)A^/2 _ 3(l + W(2))C(iVo)_ ^^^^ 

V 3ijj(2) + 5 / 3w(2) + 5 

From this we can determine 0(A^*) and hence C(-^oo) using equation (fl^. since ~ W2- Note that equation ([2^ works 
well only when the system is adiabatic. Despite this our equation for the final value of total curvature perturbation 



i.e. C(-^oo) gives accurate values as can be seen from Fig. 2.(c 
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C. Numerical results 



In order to study the applicability of these results, we have compared the analytical values with tho se ob tained 
from full numerical evolution of the equations. An example of the numerical calculations is shown in Fig. l.(a) where 
we have calculated the evolution of a system where an initially dominating matter fluid (i^(i) — 0) that decays into 
radiation (ijJ{2) — 1/3). Such a transition occurs e.g. when the inflaton decays to radiation via coherent oscillations 
or in the curvaton scenario (although it is typically assumed that the curvaton is initially subdominant, see e.g. 

0) . The initial values of the energy densities are p(i)(A^o) = 0.999, p{No) = 1.0 and the initial state is adiabatic, 
C{No) = C(i)(^o) = 1.0 with (/.(TVo) = 0.01. 

The interaction strength ratio is given by T/H^, — 100 and the transition function by f{N) = (tanh[(iV — N^,)/t] + 

1) /2. In order to model the rapid transition we choose T = 1.0 X 10~^ with - Nq = 3. 





^l+,a,n/ Cn-,num 
~4'*l Cl+,imm 



l.(a) The evolution of a matter-radiation system. 
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l.(b) Ratio of analytical to numerical values of C{l)^ 



FIG. 1: An example of a matter-radiation system with initial matter domination (a). Comparison between the numerical and 
analytical values of the jump C{i)-i- ^ function r for different values of the interaction strength, F/J/"* (b). 



Figure l.(a) shows how C(i) experiences a sudden decrease at A'^ = 3, until it slowly begins to decay. Similarly, the 
total curvature perturbation, C, begins to decline rapidly after the jump until it quickly reaches a constant value. The 
evolution of the metric perturbation, 0, is less dramatic and it reaches a constant value at a rate similar to that of 
The final state is again adiabatic, C = as expected. 

If the second fluid dominates from the beginning, like in the curvaton scenario, the behavior of perturbations is 
quite similar to the one pictured above: the jump of ^(i) depends mainly on the magnitude of T/H^ and the initial 
values. The final value of C again depends on T/H^, but it is more strongly related to the initial values of Po{i) and 
Po- 

The comparison between the analytical result, Eq. (I12p . and the numerical calculation is shown in l.(b) In the 



figure l.(b) we plot the ratio of the analytical formula to the numerical result for different values of F/iJ*. As can be 



seen from the figure, the analytical formula, Eq. (fT^ . is very good when F/i/* < 1 but it begins to over-estimate the 
magnitude for larger F/7J,. For ^/H^, > 11.75 the value of C(i)+ = — 0(A^*) is a better estimate. 

The comparisons between the full numerical calculation and the analytical approximation for asymptotic values, 
Eq. (Hn]), are presented in Fig. 2 for the same system as studied above with the first fluid initially dominant (a) 
or sub-dominant (b). In the first case the initial densities of fluids have been set to /9(i)(Afo) = 0.99, p{No) — 1.0 
and the system has been chosen to be non-adiabatic. However, since fluid one dominates the system is very close to 
the adiabatic state. In the second case where the second fluid is dominating from the beginning, we set p(i)(A'o) = 
10~^'° and p{No) = 1.0. We have calculated the evolution of the system with adiabatic and non-adiabatic initial 
perturbations. The latter of these is related to the curvaton scenario, where the initial values are (^(i)(Ao) — 1.0, 
C(2)(-^o) = 0. The metric perturbation is chosen as in the previous case. From the figure we see that the analytical 
formulae generally give a good approximation to the exact numerical calculation. Again, when T/H^, is smaller, 
the approximation given by Eqs. p2p and ([H]) is better, which is what we expect since then f{N) is more closely 
approximated by the Heaviside function. 
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2. (a) Non-adiabatic matter domination initially. 

C(JVoc,) 

0.4 



^iog(r/H,) 



Cl+, an 
,mu: 



0.5 1 1.5 2 25 3 

2.(b) Adiabatic radiation domination initially. 
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2.(c) Non-adiabatic radiation domination initially. 
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The comparison between numerical (continuous) and different analytical values (dashed) of <^{Noo), in a matter- 



FIG. 2: 

radiation system when the initial state is dominated by matter (a), adiabatic radiation (b) and non-adiabatic radiation (c) 



IV. CONCLUSIONS AND DISCUSSION 



In the present article we have discussed the evolution of perturbation variables in a two-component interacting 
fluid system. This has been studied previously by a number of authors, e.g. [1, d, 0, [13] ■ Our treatment goes beyond 
the sudden decay approximation by including all relevant perturbations explicitly and including a time scale, t, that 
describes how quickly the interaction is turned on. Such an interaction is well motivated by physical considerations 
and worth considering particularly in the fluid approach, where the microphysical properties of the fluids are otherwise 
hidden. 

We have seen how the the curvature perturbation can be efficiently suppressed in a system with interacting fluids. 
A fluid that decays rapidly, tH^ <C 1, T/H^ 3> 1, generally diminishes the curvature perturbation associated with 
that fluid which then subsequently affects the total curvature perturbation. The analysis shows that the decrease of 
perturbations depends crucially on the ratio T / in addition to the initial conditions. 

We have derived two expressions for different values of T/ which give the magnitude of this jump and compared 
these analytical results to the numerical calculations, showing good agreement. We have also derived a formula 
to approximate the final value of the total curvature perturbation, ^, after the transition. Again we find that the 
approximation is good when the decay is fast enough, i.e. T/H^ > fO~^. 

Processes with efficient suppression of the large-scale curvature perturbation can in general be realized in a scenario 
where the three physical time scales have a hierarchy, <C T, l/r. Such a hierarchy can be achieved naturally if the 
time scales are associated with different physical processes. If, on the other hand, the decay processes become effective 
only when T ~ H^, no large net effect results. Interesting avenues to explore in this respect are processes including 
phase transitions and/or different scales, making many cosmological phase transitions a potential source of interesting 
phenomena. In addition, inflationary scenarios with multiple phases can provide a suitable framework for modifying 
the large-scale perturbations, e.g. when the intermediate phases between inflationary periods are dominated by string 
networks p^ . 
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APPENDIX 



Naive integration of the equation of motion for C(i) (fTT|) over the interval [iV, — e, TV* + e], where e > 0, gives 



N,+e 



[ 9HW (C-C(i)) + -n-, {^ + C(^i))]dN. (A.i) 

By going to the limit e 0+, the term proportional to C — ^(i) clearly disappears, since the integrand is limited and 
the interval of integration goes to zero. The other term requires more care. Clearly the numerator of the integrand 
is proportional to 9'{N — N^,) (or S{N — N^,)) while the denominator has a term proportional to the 9{N — N^,). In 
order to handle this term properly, we include a C°° test function (with compact support) ip{N) into our integrals 
and handle the integrals in terms of weak convergence of distributions. 
We study integrated form of Eq. pT|) 

""'^^ C[,^iNMN)dN = [""'^^ Ep^fM}Po(^^^N) - C(i)(A^)) + f'iNmN) + C(i)(iV)))^(iV)dA^ (A.2) 



N,-<i JN. 



and after evaluations take the limit e — > 0+. We also split the (-/V)-function into two parts: the continuous part, 

C(i)(iV), and the jump AC(i)0(A^ - N,), where AQi) = C(i)(^*+) " C(i)(^*-) = C(i)+ " C(i)-- 
From l.h.s. of Eq. (jA.2|) we get 



N,+e 



h{N,)=\mi^l C[i)V{N)dN = AC(i)(p{N,). (A.3) 



and first term on r.h.s. vanish at e-limit. Next we take the the continuous part of + term on the right side of 
equation (|A.2p which yields 



= -_Bm^ (xi«.-,,».+,l^lnWW) + «(«-JV.)]') + (A.4) 

where we have introduced the function X[Ar,-e,Af,+e]j 

X[N,^e,N,+e] = 0{N -N,+e)- e{N - TV, - e), (A.5) 

and defined x{N) = 3iy(iV)(l + W(i))/r. 

The last term of (|A.4[) vanish at the limit e 0+ and we can now integrate the other term in parts when all terms 

not proportional to the derivative of X[N,-e,N,+e] gives vanishing contribution to the integral in the limit e 0+, 
too. Thus, 

hiN,) = ^{N,){\n[x4 - ln[.T, + 1])(0, + C(i),). (A.6) 
Finally, we integrate the jump-part of </i + Qx) term: 

/■^•+" , r^-+' AC(iMN)x(N)0'(N ~ N^) 

The first term in r.h.s. is clearly A(^(^ijip{N^) and the second term can be written in the form 
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Again term proportional to x'{N^) vanishes when e 0+ and the other term can be integrated by part. We obtain 

h{N^)-^ -AC(i)<^(^*)-AC(i)a;*v5(^*)(ln[x*]-ln[x* + l]). (A.8) 
Now, because Eq. (IA.2P imphes 

h=h + h, (A.9) 
by adding up all the terms and canceling all non-zero common factors we find 

^ 2 + x.ln[^] + ^^-'^^ 



where the parameter is defined by 

_ 3i/(iV,)(l + c^(i)) 



(A.ll) 
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